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Abstract 

In these proceedings, we briefly review how jets can be reconstructed in heavy-ion collisions. The main point we 
address is the subtraction of the large contamination from the underlying event background. We first present the main 
ingredients needed to define the jets and perform the background subtraction and then discuss the efficiency of the 
subtraction for different jet algorithms and background-estimation methods. 



1. Introduction 

Jets, naively seen as a proxy for the hard partons produced in a final state, have always been powerful tools in 
collider physics. Recently, interest has grown in using jets in heavy-ion (HI) collisions as a probe of the properties of 
the dense medium created in these collisions. 

The main limitation to reconstruct jets in HI collisions is the fact that they have to be isolated from the large 
Underlying Event (UE), of (9(100 - 300) GeV per unit of y - area. This UE contamination to the jet has to be 
subtracted if one wants to reconstruct its momentum correctly and obtain an unbiased measurement of the effects of 
the medium. 

In practice, both RHIC and the LHC have already accumulated a large amount of data and have tried various 
approaches to address that issue (see e.g. |lT]|2l[3]|4|). 

Here, we describe a jet-based method fSj for full jet reconstruction derived from the original proposal |6| for 
pileup subtraction in pp collisions. We shall first present the method (see eq. ([T} below) and describe its main 
elements: defining a jet, defining its area and estimating the background density per unit area. We shall then study 
using Monte-Carlo simulations what kind of performances one might expect in practice. 

Most of the effects we shall describe below can be estimated analytically. However, since this often involves rather 
technical arguments, we will not reproduce them in these proceedings and redirect the interested reader to Ref. Q. 
Furthermore, we shall concentrate on the situation at the LHC. Similar results are obtained for RHIC and can be found 
inHa. 

2. Defining background-subtracted jets 

Master formula. If the background one wants to subtract from the jets is sufficiently uniform, one can characterise it 
by two numbers: its average transverse-momentum density per unit area p and its fluctuations per unit area cr. For a 
given jet, the contamination due to the additional background will thus be pAjet + cr -y/Ajet, where Ajet is a measurement 
of the area of that jet. To correct for the background contamination, we can thus define a (4-vector) subtracted jet 
momentum as 

(1) 




ELSEVIER 
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In order to apply that master formula, three building blocks, described below, are necessary: define jets, define their 
area and estimate the background density p. 

Before getting to that, let us mention that eq. ([T} is incomplete in two major aspects. The first one comes from 
the fluctuations of the background, leaving an uncertainty of the order of cr -y/Aj^. The second source is referred to 
as back reaction fT\: since background particles have a non- vanishing momentum, they can affect the clustering of 
the "hard" particles; this means, on top of the pure background contamination (pAjet), the hard contents of the jet can 
change when considering it together with the background. We will show some examples of these effects in Section[3] 

Jet definitions. Defining jets from a list of input objects — e.g. particles or calorimeter towers — has been the 
topic of many discussions over the last few decades (see e.g. |8j). Here, we shall concentrate on recombination-type 
algorithms that define the distances 

dij = min(^?J,^,2p(A4 + Ay2.), (2) 
diB = kl^R^ (3) 

respectively between any pair of particles and for a single particle. The clustering works by successively identifying 
the smallest distance; if it is a djj recombine particles / and j, otherwise, call / a jet. 

The geometric factor in dij basically ensures that collinear particles will be clustered in the same jet. The k,- 
dependent prefactor varies from one algorithm to another: the k, algorithm ||9] [TOl [1 1 1 has p = 1 so that soft QCD 
emissions are clustered early in the sequence; the simple case p = corresponds to the Cambridge/Aachen (C/A) 
algorithm IIT2l[T3l . useful e.g. for jet substructure studies; the anti-A:, algorithm |14 |, p = -1, will cluster jets around 
hard objects, producing hard jets with rigid, circular, shapes. Note that the anti-A:, algorithm is the default choice for 
most of the LHC experiments. All these algorithms are accessible using the FastJet package ifTSl fT&l and we shall fix 
R = 0.4 in what follows. 

Recently, a lot of progress has been made that involve making use of the substructure of the jets. To illustrate 
that, we shall consider, on top of the 3 algorithms introduced above, the C/A algorithm supplemented with a filter 
(C/A(filt)) IfTTl . The latter works by reclustering each individual jet with a smaller radius (we shall use R/2) and only 
keep the 2 hardest subjets. Because of the collinear nature of QCD branchings, the filtering procedure is expected to 
keep most of the QCD part of the jet and reject part of the Underlying-Event contamination. 

Jet areas. Practically, the area of a jet is defined Q as the region in which it captures (infinitely) soft particles (ghosts). 
Here, we will use the active area of jets determined by adding a dense coverage of ghosts to the evenj^ each carrying 
a "quantum" of area, clustering them together with the jets, and obtain the area by summing the ghost it contains. This 
procedures mimics the addition of a background to the event with the exception that ghosts are infinitely soft. 

The precise details of how the jet area is defined are largely irrelevant. However, we would like to stress that for 
the computation of the background density (see below), we advise the use of the explicit_ghost option of FastJet ifTSl 
to correctly handle empty regions of the event. 

Background estimation. In f6l|, the following method has been suggested: first cluster the event with the kt or C/A 
algorithm with a radiua^/?p to use, as an estimator for the background density per unit area. 



p - median 




where the computation of the median includes all the jets up to a maximal rapidity jmax- We refer to that choice as the 
global range below. 

To better take into account the non-uniformities of the background, we have realised that there is some interest in 
using only jets in the vicinity of the jet we want to subtract: 

Pm = median ] [ , (5) 



This will not modify the clustering of the hard jets, provided the algorithm is infrared-safe. 

^This choice of jet definition for the estimation of the background is independent of the choice made to reconstruct the jets in the event. For the 
former, we suggest using the kt or C/A algoiithm as it avoids having many jets with a small area. 
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Figure 1 : Distribution of Ap, (red liistograms) at RHIC for each of our 4 jet algorithms, together with a Gaussian (black curve) whose mean (solid 
vertical line) and dispersion are equal to (Ap,) and cTAp, respectively. 



where the subscript 'R(\) explicitly refers to the fact that we use a local range centred on the jet j for which we want 
to estimate the backgrouncj^ 

We have considered 3 different options: a CircularRange{A) keeping jets within a distance A (we used A - 3R) of 
the jet, a DoughnutRangeiS, A) keeping jets with a distance from the jet between 6 and A (we used S - R, and A -2R 
or 3/?) and a StripRange(A) being a rapidity strip extending between yjet - A and yjet + A (we used A = 2/? or 3R). 

The use of a median in Q and (j5]l rather than an average is made to reduce the effect of the hard jets on the 
estimation of p. In order to further reduce that bias, we have also tested the exclusion of the two hardest jets in the 
event from the jets used for the median computation. 



3. Expected performances 

Framework. In this Section, we study how the method presented above performs in practice. 

To achieve that, we shall embed pp hard events in HI background, apply our background-subtraction method and 
see how precisely the original hard jets are reconstructed. The pp events will be generated using Pythia 6.4 |18| and 
for the AA events we have used Hydjet 1.6 |fT9l . We note that, for 0-10% central collisions, the HI background has an 
average density (p> of about 100 GeV at RHIC and 310 GeV at the LHC, with fluctuations <cr) around 8 and 20 GeV, 
respectively. 

Jets are reconstructed and subtracted for both the hard event alone and the hard event embedded in the heavy-ion 
background (which we shall refer to as the full event). We have concentrated on the two hardest jets in the hard event. 



^Note that if the use of a local range allows to better estimate of the background at the position of the jet we want to subtract, it may also have 
the inconvenient to increase the sensitivity to quenching effects. 
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Figure 2: Effect of the choice of range on the average p, shift, Ap,, as defined in eq. J6j, for LHC kinematics. In this figure and those that follow, 
the yellow band con'esponds to 1% of the p, of the hard jet. The label "2 excl" means that the two hardest jets of the event have been excluded 
from the median computation. 



For each of them, we search for a matching fulljet in the full event by requiring that the hard contents of the full jet 
accounts for at least 50% of the original hard jejj 

For a matching pair of jets, we consider the difference between the subtracted fulljet and the subtracted hard jet: 

A AA,sub pp,sub 

^Pt = P, - Pt ■ (6) 

An efficient background subtraction would directly translate into small Ap,. Though we could obtain the full 
distribution of A/?,, see Fig. [T]for RHIC kinematics, they can be considered as sufficiently Gaussian for the purpose 
of comparing different background-estimation ranges and jet-reconstruction algorithms, and we shall concentrate on 
its average (Ap,) and dispersion 

An effective subtraction would then be characterised by small values of (A/?,) and cr^p,. 

For the results presented below, the k, algorithm with a radius Rp - 0.5 has been used for the median-based 
estimation of the background density. 

Choice of range. We start our investigation by looking at how the choice of the range made to estimate the background 
density affect the subtraction. On Fig. [2j we see that for all the ranges we consider, the average uncertainty on the 
reconstructed jet transverse momentum is < 1%, with a small preferenc^for local ranges at large p,. 

In choosing the range, there is a tension between remaining in the vicinity of the jet one wants to subtract and 
having enough jets in the range to estimate the median properly. As a rule of thumb, it can be shown that a local range 
needs to contain at least 9 jets to give a reasonable estimate of p. 

For what follows, we shall use the DoughnutRange(/?, 3/?) as a default choice. 

Choice of algorithm. We now focus on the dependence on the specific choice of algorithm and first consider the 
average shift (Ap,). It is plotted on Fig. [Sjleft) for the four algorithms under consideration. We see that while the 
anti-^, and C/A(filt) algorithms give an average shift close to zero, the k, and C/A algorithms have a larger, negative, 
shift. 

Since the estimation of the background density p is common to all 4 algorithms, this can actually be traced back to 
a difference in back-reaction as seen on Fig. [3](right). For all except C/A(filt), the differences observed in the average 



'^Though we do not show it here, the matching efficiencies depend slightly on the jet definition and are typically above 95% for /), > 10 GeV at 
RHIC and above 98% for p, > 40 GeV at the LHC. 

'For a limited rapidity coverage, a global range behaves as a local range. However, at large p,, jets tends to be produced at smaller rapidity 
where the background is larger and a global estimation of p leads to an under-subtraction. 
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Figure 3: Left: average sliift (Ap,), as a function of Pf.hardi shown for the LHC. Right: contribution to (Ap,> due to back-reaction. Note that the 
back-reaction results ai'e independent of the range used for estimating p in the heavy-ion event. 
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Figure 4: Dispersion (T^pj for the LHC kinematics. 



shift correspond to differences in back-reaction and the fact that the anti-A;, has a small shift is a direct consequence of 
its rigidity. In the case of the C/A(filt) the final small shift can be explained by a (fortuitous) cancellation between the 
offset due to back-reaction and a filtering bias. 

Finally, we turn to the dispersion cr^p,. The observed dispersions for the LHC are shown on Fig. |4] While anti-fc, 
was the more robust in terms of average shift, when we consider the dispersion, C/A(filt) shows better performances. 
What we really see here is precisely the effect of filtering which practically reduces the area of the jet, thus reducing 



the dispersion which is proportional to 



With our choice of parameters, the area after filtering is about 1/2 of 



the original C/A-jet area and we may thus expect a decrease of a factor 1 / V2 in the dispersion, which is consistent 
with what is observed on the figure at low p,. 

Additional remarks. To conclude this Section, we want to make a few remarks on other tests that can be performed 
(see for more details). 

First, to test the robustness of our results, we have checked the effect of quenching on the embedded hard events. 
This can been done e.g. by running PyQuen |20| on the Pythia jets. We have noticed a negligible effect on the results 
obtained with the anti-A:, algorithm and an effect smaller than 2% for the C/A(filt) case. 

Then, we can investigate if our conclusions hold for more peripheral collisions, where flow effects may increase 
the fluctuations. For the average quantities presented here, we have checked that our conclusions remain unchanged 
in the case of non-central collisions. However, it is interesting to mention that if we study the average shift (Ap,) as a 
function of the azimuthal angle between the jet and the reaction plane, we see oscillations of a few GeV. The use of a 
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range with a finite azimuthal coverage, like the DoughnutRange, reduces these oscillations compared to a StripRange, 
but further developments would be needed to reduce them completely. 

4. Conclusions 

To conclude, we have studied the possibility to use the jet-area-based subtraction method initially proposed in ||6l 
to reconstruct full jets in heavy-ion collisions. 

We have found that the use of local ranges to reduce the non-uniformities of the background (like its rapidity 
dependence) improves the background subtraction. It would be interesting to see if using different ranges (bearing in 
mind that they should contain at least ~9 jets to give a reliable estimate of p) could help obtaining an unceratinty on 
the estimation of p. 

We have studied subtraction performances by embedding Pythia events in Hydjet heavy-ion events and measured 
for the two hardest jets in each event Ap,, the transverse momentum difference between the original (subtracted) hard 
jet and the embedded/full jet after subtraction. While knowledge of the full Ap, distribution is relevant, e.g. for 
unfolding purposes 121], we have chosen to focus only on its average and dispersion. 

For the average shift {Ap,}, we have observed that the anti-A:, (because of its rigidity) and the C/A algorithm 
supplemented with a filter give a result close to zero, though one has to notice that for the latter, this is the result of a 
cancellation between two effects. 

In terms of dispersion the C/A(filt) algorithm, which has a smaller sensitivity to the UE, as a consequence of its 
reduced area, shows better performances than the other algorithms. 

Finally, we have noticed (see |5 | for more details) that these results also hold when quenched jets are embedded 
or when one considers more peripheral collisions. In this last case, there is a residual shift when one considers the 
shift as a function of the azimuthal angle between the jet and the reaction plane that is left for further studies. We also 
note that, though the C/A filt algorithm has a smaller dispersion, it may be more affected by quenching than the anti-A;, 
algorithm (e.g. because of jet broadening) and it may thus be helpful to consider both to obtain the most complete and 
valuable information. 
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